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Abstract 

A  new  approach  for  design,  generation,  computerized  simulation  of  meshing  and  contact,  and  stress  analysis  of  fonmte 
cut  spiral  bevel  gear  drives  is  reported.  The  approach  proposed  is  based  on  application  of  four  procedures  that  permit  in 
sequence  to  provide  a  longitudinally  directed  bearing  contact,  a  predesigned  parabolic  function  of  transmission  errors,  a 
limit  to  the  shift  of  bearing  contact  caused  by  errors  of  alignment,  and  to  automatically  generate  the  fimte  element  models 
for  application  of  a  general  purpose  finite  element  analysis  computer  program.  A  new  computer  program  has  been 
developed  that  provides  optimized  machine-tool  settings  for  finishing  of  spiral  bevel  gears  utilizing  the  Tormate 
grinding  process  with  minimized  transmission  errors.  The  new  computer  program  is  based  on  the  four  procedures 
mentioned  above.  The  developed  theory  is  illustrated  with  several  examples  of  design  and  computations. 
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Nomenclature 

Profile  angle  of  parabolic  blade  for  gear  head  cutter  at  mean  point  M 
Blade  angle  for  straight  blade  of  pinion  head  cutter 
Elastic  approach  of  contacting  surfaces 

Pitch  cone  angles  of  pinion  (/  =  1)  and  gear  (/  =  2) ,  respectively 
Machine  root  angles  for  pinion  (i  =  1)  and  gear  (/  =  2) ,  respectively 
Root  cone  angles  for  pinion  (i  =  1)  and  gear  (/  =  2) ,  respectively 
Shaft  angle 

Tangent  to  the  path  of  contact  on  the  pinion  (/  =  1)  and  gear  (/  =  2)  surface,  respectively 
Surface  parameters  of  the  pinion  and  gear  head  cutters,  respectively 
Surface  parameters  of  the  pinion  and  gear  fillet  parts  of  the  head  cutter 
Surface  parameters  of  the  pinion  and  gear  fillet  parts  of  the  head  cutter 
Angle  formed  between  principal  directions  (i,y  =  1, 2,p,g) 

Angle  of  rotation  of  the  pinion  (i  =  1)  or  gear  (/  =  2)  in  the  process  of  meshing 

Angle  of  rotation  of  the  cradle  in  the  process  of  generation  of  the  pinion 

Angle  of  rotation  of  the  pinion  in  the  process  of  generation 

Fillet  radii  for  the  pinion  and  the  gear 

Angular  velocity  of  the  pinion  (in  meshing  and  generation) 
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Pinion  (/  =  1)  or  gear  (/  =  2)  tooth  surface 

Pinion  {k  =  p)  or  gear  {k  =  g)  generating  surface 

Function  of  transmission  errors 

Pinion  and  gear  axial  displacements,  respectively 

Errors  of  the  offset  and  shaft  angle,  respectively 

Parabolic  coefficient  of  gear  head  cutter  blade 

Dedendum  of  pinion 

Mean  cone  distance  and  face  width 

Coefficients  of  modified  roll  for  pinion  generation  (/  =  1, 2, 3) 

Unit  vectors  of  principal  directions  of  pinion  and  gear  tooth  surfaces,  respectively 
Unit  vectors  of  principal  directions  of  the  cutter  surface  for  pinion  and  gear,  respectively 
Blank  offset  for  pinion 
Normal  curvature  of  surface  Zy 
Relative  normal  curvature 

Principal  curvatures  of  pinion  and  gear  tooth  surfaces,  respectively 
Principal  curvatures  of  cutter  surfaces  for  pinion  and  gear,  respectively 
Second  derivative  of  transmission  function  ^2(^) 

Gear  ratio 
Mean  contact  point 

Matrices  of  coordinate  transformation  from  system  5,  to  system  Sj 
Unit  normal  and  normal  to  surface  Z^  represented  in  coordinate  system 
Tooth  number  of  pinion  (/  =  1)  or  gear  (/  =  2) 

Installment  angle  for  head  cutter  of  the  pinion 

Radii  of  the  head  cutter  at  mean  contact  point  for  the  pinion  and  gear 

Gear  cutter  radius 

Position  vector  represented  in  system  S.  (z  = 

Surface  parameter  of  the  pinion 

Coordinate  system  (z  =  1, 2,fcj,Z>2,^,/,/W|,W2,/7,g) 

Radial  setting  of  the  head-cutter  of  the  pinion 
Gear  tooth  surface  parameter 

Velocity  of  contact  point  in  its  motion  over  surface  Z^ 

Transfer  velocity  of  contact  point  in  its  motion  with  surface  Z^ 

Relative  velocity  at  contact  point  (z,y  =  l,2,Cj ,/?) 

Components  of  the  velocity  of  the  contact  point  in  its  motion  over  Z^ 

Sliding  base  for  pinion 
Machine  center  to  back  for  pinion 

Center  distances  of  fillet  circular  arcs  of  pinion  and  gear,  respectively 
Machine  center  to  back  for  gear 
Length  of  major  axis  of  contact  ellipse 
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Chapter  1 

Developed  Theory 


1.1  Introduction 

Design  and  generation  of  spiral  bevel  gears  is  a  subject  of  intensive  research  of  scientists,  designers,  and 
manufacturers.  High  quality  equipment  such  as  CNC  (computer  numerically  controlled)  machines  and  tools  have  been 
produced  by  Gleason  Works  (USA)  and  Klingelnberg-Oerlikon  (Germany  and  Switzerl^d)  for  precise  manufacturing. 
However,  the  quality  of  spiral  bevel  gears  depends  substantially  on  the  design  of  machine-tool  settings  that  have  to  be 
determined  for  each  design  of  spiral  bevel  gears.  The  machine-tool  settings  are  not  standardized  and  their 
determination  is  a  subject  of  high  quality  research.  The  freedom  of  choosing  optimized  machine-tool  settings  allows  us 
to  generate  low-noise  stable  bearing  contact  spiral  bevel  gears. 

The  subject  of  this  research  project  is  the  enhanced  design  of  formate  cut  spiral  bevel  gears.  Formate  is  a  Gleason 
Works  trademark  applied  for  nongenerated  spiral  bevel  gears.  The  term  “formate  cut”  means  that  the  gear  is  held  at 
rest  and  grinding  or  cutting  wheel  (head  cutter)  does  not  orbit  the  cradle.  Thus,  the  gear  tooth  surface  is  a  copy  of  the 
tool  surface  that  is  a  surface  of  revolution  (see  Section  1.3).  The  pinion  tooth  surface  is  generated  by  a  head-cutter 
while  the  tool  and  the  pinion  perform  related  rotations.  The  pinion  tooth  surface  is  generated  as  the  envelope  to  the 
family  of  the  tool  (head-cutter)  surfaces.  The  advantage  of  the  formate-cut  method  is  the  high  productivity  of  gear 
generation.  The  main  problem  of  the  design  of  formate-cut  spiral  bevel  gears  is  the  conjugation  of  tooth  surfaces  of  the 
gear  and  the  pinion. 

The  contents  of  the  report  cover:  (i)  computational  procedures  of  the  proposed  approach  that  permit  the  design  of  a 
longitudinal  bearing  contact  and  a  parabolic  function  of  transmission  errors,  and  investigation  of  the  influence  of  errors 
of  alignment;  (ii)  algorithms  of  computer  programs  developed  for  synthesis,  analysis  of  meshing  and  contact,  and 
automatic  development  of  finite  element  models;  (iii)  results  of  computations  for  several  cases  of  design  of  formate  cut 
spiral  bevel  gear  drives;  and  (iv)  results  of  finite  element  analysis  of  several  design  cases  of  both  face-milled  generated 
and  formate  cut  spiral  bevel  gear  drives. 

1.2  Basic  Ideas  of  Proposed  Approach 

Localization  of  Contact  The  contact  of  pinion-gear  tooth  surfaces  is  localized  due  to  the  mismatch  of  generating 
surfaces  (surfaces  of  the  head-cutters)  that  is  achieved  by  application  of  different  dimensions  and  profiles  of  the  gear- 
pinion  head-cutters. 

The  generated  pinion-gear  tooth  surfaces  are  in  point  contact  at  every  instant.  The  tooth  contact  under  load  is  spread 
over  an  elliptical  area  and  the  bearing  contact  is  formed  as  a  set  of  contact  ellipses.  The  magnitude  and  orientation  of 
instantaneous  contact  ellipse  requires  information  about  the  elastic  deformation  of  tooth  surfaces  and  the  principal 
curvatures  and  the  principal  directions  of  tooth  surfaces  at  each  instantaneous  point  of  their  tangency. 

It  will  be  shown  later  (Section  1.5)  that  the  path  of  contact  on  gear  tooth  surface  is  oriented  longitudinally.  Such  a  path 
of  contact  is  .favorable  for  increasing  the  contact  ratio  (of  the  average  number  of  teeth  that  are  in  contact 
simultaneously),  allows  edge  contact  to  be  avoided  and  improves  the  conditions  of  transfer  of  meshing  wherein  one 
pair  of  teeth  is  changed  for  the  neighboring  one. 
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The  TCA  (Tooth  Contact  Analysis)  computer  program  developed  by  the  authors  allows  simulation  of  meshing  of 
surfaces  and  investigation  of  the  shift  of  bearing  contact  caused  by  errors  of  alignment.  In  some  cases  of  design,  it  may 
become  necessary  to  adjust  the  bearing  contact  -  to  deviate  it  from  the  longitudinal  direction  -  to  reduce  the  shift  of 
bearing  contact  caused  by  misalignments. 

Application  of  Parabolic  Profiles  of  Blades  of  Gear  Head-Cutter.  We  have  mentioned  above  that  the  gear  tooth 
surface  is  generated  as  the  copy  of  the  surface  of  the  gear  head-cutter.  Usually,  straight  profile  blades  are  applied  for 
the  gear  head-cutter.  Such  blades  form  the  head-cutter  generating  surface  as  a  cone. 

Application  of  parabolic  profiles  of  the  blades  instead  of  straight  line  profiles  enables  to  modify  the  gear  tooth  surface 
in  favor  of  conjugation  of  pinion-gear  tooth  surfaces. 

Simulation  of  Meshing.  The  algorithm  of  computerized  simulation  of  meshing  provides  continuous  tangency  of 
contacting  gear  tooth  surfaces  Z|  and  E2  as  the  equality  of  position  vectors  and  unit  normals  nj^*^  =  n^^^ 

(see  details  in  Section  1.8).  Instead  of  equality  of  unit  normals  to  contacting  surfaces,  collinearity  of  surface  normals 
may  be  considered.  Position  vectors  r^*^(/  =  l,2)  and  unit  normals  nj,'^(/  =  l,2)  are  considered  in  fixed  coordinate 
system  (Fig.  1.1). 


Figure  1.1:  Illustration  of  tangency  of  surfaces  Ej  and  E2  . 

Application  of  TCA  permits  the  determination  of  the  paths  of  contact  on  pinion-gear  tooth  surfaces,  dimensions  and 
orientation  of  instantaneous  contact  ellipse  at  each  point  of  contact,  the  transmission  function  ^(^)  that  relates  the 
angles  of  rotation  of  the  gear  and  the  pinion,  and  function  A^(^)  of  transmission  errors  that  is  determined  as 

=  (1.1) 

J\2 

where  and  N2  are  the  numbers  of  teeth  of  the  pinion  and  the  gear. 

Predesign  of  Parabolic  Function  of  Transmission  Errors.  The  output  of  TCA  for  a  misaligned  gear  drive  shows  that 
the  transmission  function  ^{0^)  is  a  piecewise  discontinuous  almost  linear  function  (Fig.  1.2(a)).  The  transfer  of 
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meshing  with  the  transmission  function  (Fig-  1.2(a))  and  function  A^((^[)  of  transmission  errors  (Fig.  1.2(b)) 

is  accompanied  with  the  jump  of  the  angular  velocity  a>2  and  therefore  noise  and  vibration  are  inevitable. 


(a) 


I  M  I 

(b) 

Figure  1.2:  Transmission  function  (a)  and  function  of  transmission  errors  (b)  of  a  misaligned  gear  drive. 

The  purpose  of  synthesis  of  low-noise  gear  drives  is  to  transform  the  transmission  function  ^(^)  as  shown  in  Fig. 
1.3(a)  and  obtain  a  parabolic  function  of  transmission  errors  shown  in  Fig.  1.3(b).  The  modified  function  A^(^)  (Fig. 
1.3(b))  is  a  continuous  one  and  is  determined  as 

where  a  is  the  parabola  coefficient. 

It  is  important  to  recognize  that  the  flmction  of  transmission  errors  must  be  negative  and  therefore  the  gear  will  lag 
with  respect  to  the  pinion.  Therefore,  the  transfer  of  meshing  is  accompanied  with  elastic  deformations,  and  as  a  result, 
the  gear  ratio  is  increased. 

The  desired  parabolic  function  of  transmission  errors  is  provided  by  the  respective  synthesis  of  the  gear  drive.  The  idea 
of  such  a  synthesis  is  based  on  the  predesign  of  a  parabolic  function  obtained  by  application  of  modified  roll  during  the 
process  of  generation.  Modified  roll  is  provided  due  to  application  of  nonlinear  relations  between  the  angles  of  rotation 
of  the  cradle  of  the  generation  machine  and  the  pinion  during  the  process  of  pinion  generation  (see  Section  1 .5). 

The  advantage  of  the  predesign  of  a  parabolic  function  of  transmission  errors  is  the  absorption  by  such  a  function  of 
the  linear  discontinuous  function  of  transmission  errors  caused  by  misalignment.  This  is  illustrated  by  drawings  of  Fig. 
1.4(a)  that  shows  a  linear  function  A^P^(^)  =  b^  (Fig-  l-4(n)).  Parabolic  function  A^2^(^)  is  predesigned  and 
linear  function  A^P^(^)  is  caused  by  misalignment.  The  sum  of  functions  and  A^P^(^)  is  obtained  as  a 

parabolic  function  and  represented  in  Fig.  1.4(a)  as  function  A^2(l^i)  “  •  The  parabola  coefficient  a  is  the  same 

in  functions  and  At//2iWi)- 
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(a) 


(b) 


Figure  1.3:  Transmission  function  (a)  and  function  of  transmission  errors 
(b)  of  a  gear  drive  with  a  predesigned  parabolic  function. 


The  drawings  confirm:  (i)  the  predesigned  parabolic  function  absorbs  indeed  linear  function  ,  and  (ii)  the 

resulting  function  Aw2('^i)  is  slightly  dislocated  with  respect  to  origin  of  ,  function  A^2(V'i)  is  an 

asymmetrical  one.  Parameters  c  and  d  of  the  origin  of  Aif/^iyfO  are  determined  as  d  =  b^  lAa,  c  =  A/2a.  In  the 
process  of  meshing,  wherein  several  cycles  of  meshing  are  considered  (Fig.  1.4(b)),  the  resulting  function  Ay/2{W\)  of 
transmission  errors  becomes  a  symmetric  parabolic  function. 

The  authors  have  developed  computerized  synthesis  of  formate-cut  spiral  bevel  gears  that  is  based  on  simultaneous 
application  of  local  synthesis  and  simulation  of  meshing  and  contact  of  gear  tooth  surfaces. 

The  purpose  of  local  synthesis  is  determination  of  the  pinion  machine-tool  settings  considering  as  given:  (i)  the  gear 
machine-tool  settings,  and  (ii)  the  conditions  of  meshing  and  contact  at  the  mean  contact  point  M  such  as  the  tangent 
to  the  path  of  contact  at  M ,  the  major  axis  of  the  contact  ellipse  at  Af ,  and  the  derivative  of  the  function  of 
transmission  errors. 

The  purpose  of  simulation  of  meshing  and  contact  is  the  determination  of  the  bearing  contact  and  function  of 
transmission  errors  knowing  the  pinion  and  gear  machine-tool  settings. 

The  applied  computational  procedures  are  iterative  processes  (see  Section  1.5)  with  the  goals  to  obtain:  (i)  a 
longitudinally  directed  bearing  contact,  (ii)  a  parabolic  function  of  transmission  errors  with  limited  value  of  maximal 
errors,  and  (iii)  reduced  sensitivity  of  the  gear  drive  to  errors  of  alignment. 
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Figure  1.4:  Interaction  of  predesigned  parabolic  function 

linear  function  caused  by  misalignment. 


1.3  Derivation  of  Gear  Tooth  Surface 

Applied  Coordinate  Systems.  Coordinate  system  is  the  fixed  one  and  it  is  rigidly  connected  to  the  cutting 

machine  (Fig.  1.5).  Figure  1.5(a)  and  (b)  show  two  sets  of  coordinate  systems  applied  for  generation  of  left-  and  right- 
hand  gears,  respectively.  Coordinate  system  S2  is  rigidly  coimected  to  the  gear.  Coordinate  system  5'g  is  rigidly 
connected  to  the  gear  head-cutter.  It  is  considered  that  the  gear  head-cutter  is  a  surface  of  revolution,  and  the  rotation 
of  the  head  cutter  about  the  axis  does  not  affect  the  process  of  generation. 

The  installment  of  the  tool  on  the  machine  is  determined  by  parameters  H2  and  Fj,  that  are  called  horizontal  and 
vertical  settings.  Parameters  Xq  and  represent  the  settings  of  the  gear. 


NASA/CR— 2003-2 12336 


7 


(a) 


(b) 

Figure  1.5:  Machine  tool  settings  for  gear  generation:  (a)  for  left-hand  gear  and  (b)  for  right-hand  gear. 

Head  Cutter  Surface.  Blades  of  parabolic  profile  (Fig.  1.6(a))  are  applied  for  the  gear  head-cutter.  Each  side  of  the 
blade  generates  two  sub-surfaces.  The  segment  a  of  the  parabolic  profile  generates  the  working  part  of  the  gear  tooth 
surface.  The  circular  arc  of  radius  generates  the  fillet  of  the  gear  tooth  surface.  The  generating  surfaces  of  the  head 
cutter  are  formed  by  rotation  of  the  blade  about  the  Xg  axis  of  the  head-cutter  (Figs.  1 .6(b)  and  (c));  the  rotation  angle 
is  0g.  Therefore,  the  generating  surfaces  are:  (i)  surface  of  revolution  formed  by  rotation  of  the  blade  of  parabolic 
profile  (part  a  ),  and  (ii)  the  surface  of  the  torus  formed  by  rotation  of  the  circular  arc  profiles  (part  b).  A  point  on  the 
generating  surface  is  determined  by  parameters  Ug  and  dg  for  the  working  surface,  and  by  and  0g  for  the  fillet 

surface.  Angle  Og  is  formed  between  the  tangent  line  of  the  blade  at  point  M  and  the  vertical  centerline  of  the  blade. 
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Parameter  measured  from  point  M  in  the  chosen  direction  is  considered  as  a  positive  one,  and  angles  and 
as  the  acute  ones.  The  apex  of  the  parabola  is  located  at  point  M  determined  by  parameter  Sg^ . 


Figure  1.6:  Illustration  of  gear  cutter  blades  and  generating  surfaces. 


The  surface  of  revolution  and  the  torus  surface  of  the  head-cutter  are  designated  as  parts  a  and  b  of  the  head-cutter 
generating  surface.  Surface  2^"^  of  the  head-cutter  is  represented  by  vector  function  t'^^{Ug,6g)  as  follows 


-  Mg  cosag  +  sinttg  -  Sg^  cosag 

{Rg  ± Sg^  sinag  +  Ug  sinag  ± a^u^g  cosag) sin ^g 

(/?g  ±  5g^  sinag  ±  Ug  sinag  ±  a^Ug  cosag)cos0g 


(1.3) 


where  Ug  and  0g  are  the  surface  coordinates,  ag  is  the  blade  angle  at  point  M ,  Rg  is  the  cutter  point  radius  (Fig. 
1.6)  and  given  by 


(1.4) 


The  upper  and  lower  signs  in  equations  (1.3)  and  (1.4)  correspond  to  the  concave  and  convex  sides  of  the  gear  tooth. 


respectively. 


The  unit  normal  to  the  gear  generating  surface  is  represented  by  the  equations 
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Equations  (1.3)  and  ( 1 .5)  yield 


*  *  *  *  dUg  80^ 


±  sina^  ±  2a cosa^ 


n'  ^(w  )=  (cosor  -2a^M  sina  )sin0  UJl^Aa^u 


^llMlA.g  /ifglMVTg 

(cos  Og  -2a^Ug  sin  Og )  cos  6g 


Surface  is  represented  in  S  as 


where 


-p^(l-cosA„,) 

rf''(X^.,eg )  =  (Xg±p^  sin  A„,)sin  0g 
(X  ±  p^.s\nX„)cos6 


i^g-^g  +p„,(l-sinag)/cosag 


Here,  p^  is  the  fillet  radius  of  the  gear  generating  surface. 

The  unit  normal  to  the  gear  generating  surface  is  represented  by  the  equations 


nW(4,^^)  = 


_nw 

ImWI’ 


Equations  ( 1 .8)  yield 


arf) 

*  5A„,  80^ 


±cosA. 


Vk^P(X„,0g)=^  sinA„,sin^g  (1.9) 

sin  QosOg 

Equations  of  Head-Cutter  Surfaces  in  ^2 .  Applying  coordinate  transformation  from  Sg  to  S2 ,  the  gear  tooth 
surfaces  can  be  represented  in  ^2  by  the  equations 

4‘’\ug,0g)  =  M,gr^‘’\ug,0g)  (1.10) 

r^'’\K>0g)  =  M2gr^g‘‘\A^,0g)  (1.11) 


M,„=1V! 


'10  0  0 
0  1  0  ±^2 
0  0  1  //j 

0  0  0  1 
rcosy„,  0  - 


M2„,= 


cosy„,  0  -siny„^ 
0  1  0 

sinv  0  cosy 


Here  V2,  Hj,  and  X(;(Fig.  1.5)  are  the  gear  machine  tool  settings.  The  upper  and  lower  signs  in  fi-ont  of  V2 
correspond  to  right-hand  and  left-hand  gears,  respectively.  The  whole  set  of  gear  machine  tool  settings  is  represented 
in  Table  1.  We  remind  that  the  generated  gear  tooth  surface  is  a  copy  of  the  surface  of  the  head-cutter,  that  is  a  surface 
of  revolution.  The  rotation  of  the  head-cutter  about  the  Xg  axis  is  necessary  for  the  cutting  or  grinding  process  but 
does  not  affect  the  shape  of  gear  tooth  surfaces. 
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_ Table  1.1:  Gear  machine  tool  settings. _ 

Blade  angle  (Fig.  1.6) 

Blade  parabolic  coefficient  _ _ _ 

Parabola  apex  location 

Cutter  (grinding  wheel)  radius  7?^  (mm)  (Fig.  1.6) _ 

Point  width  (mm)  (Fig.  1.6) 

Cutter  point  radius  (mm)  (Fig.  1 .6)  {R^-R^±P^^I2) 

Horizontal  setting  H2  (mm)  (Fig.  1.5) _ _ 

Vertical  setting  V2  (mm)  (Fig.  1.5) _ _ 

Machine  center  to  back  Xq  (mm)  (Fig.  1.5) 

Machine  root  angle  (Fig,  1.5) 

Fillet  radius  (mm)  (Fig.  1.6)  _ 


A  3D  geometric  model  of  the  gear  tooth  surfaces  is  shown  in  Fig.  1 .7. 


Figure  1.7:  3D  geometric  model  of  gear  tooth  surfaces. 


1.4  Derivation  of  Pinion  Tooth  Surface 

Applied  Coordinate  Systems.  Coordinate  systems  applied  for  generation  of  pinion  are  shown  in  Figs.  1.8  and  1.9. 
Figures  1.9(a)  and  (b)  correspond  to  generation  of  right-  and  left-hand  pinions,  respectively.  Coordinate  systems  , 

are  the  fixed  ones  and  they  are  rigidly  connected  to  the  cutting  machine.  The  movable  coordinate  systems  5j 
and  are  rigidly  connected  to  the  pinion  and  the  cradle,  respectively.  The  coordinate  systems  and  are  rotated 
about  the  axis  and  axis,  respectively,  and  the  angles  of  rotation  are  related  by  a  polynomial  function  y/x{y/c^ ) 
if  modified  roll  is  applied  (see  below).  The  ratio  of  instantaneous  angular  velocities  of  the  pinion  and  the  cradle  is 
defined  as  mxXy/x(ii/c,))  =  ■  The  magnitude  mx^{y/{)  at  is  called  ratio  of  roll.  Parameters  > 

,  and  are  the  basic  machine  tool  settings  for  pinion  generation. 
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Head  Cutter  Surfaces.  The  pinion  generating  surfaees  are  formed  by  surface  and  generated  by  straight-line 
and  circular  arc  parts  of  the  blades  (Fig.  1.10).  Surface  Z^^  is  represented  as 

(RpTSpSmap)cosdp 

r^^\sp,0p)=  {Rp+Spsmap)sindp  (1-13) 

-SpCosap 

where  Sp  and  0p  are  the  surface  coordinates;  Op  is  the  blade  angle;  Rp  is  the  cutter  point  radius  (Fig.  1.10).  The 
upper  and  lower  signs  in  equations  (1.13)  correspond  to  the  convex  and  concave  sides  of  the  pinion  tooth  respectively. 
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Equations  (1.13)  and  (1.14)  yield 


cosa^  cosOp 
cosap  sin  Op 
+  sina^ 


(1.15) 


For  the  fillet  surface  ,  we  obtain 


where 


(.Yy-  +/7ySinAy)COS0^ 
(A'y  +/7y  sin/ly)sin^^ 
--Py(l“COSAy) 


A"y  =  ±py(l-“Sinay)/cosay 

and  py  is  the  radius  of  the  tool  fillet. 


(1.16) 


The  unit  normal  to  the  pinion  generating  surface  of  part  b  is  represented  by  the  equations 


KT 


= 


ari*> 


d0p 


Equations  (1.17)  yield 


sin  Ay  cos 
sin  Ay  sin 

+  COSAy 


(1.17) 


(1.18) 


The  head-cutter  is  mounted  on  coordinate  system  (called  the  cradle  of  the  cutting  machine)  and  its  installments  are 
determined  by  settings  5^.  and  (Fig.  1.9).  In  the  process  for  pinion  generation,  coordinate  systems  5^^  and 
perform  related  rotations  with  respect  to  .  Then,  a  family  of  pinion  head-cutter  surfaces  is  generated  in  coordinate 
system  5j  determined  as 


Here  is  the  generalized  parameter  of  motion,  and  the  used  matrices  are  determined  as 


iP 


M 


ciP 


10  0  cos^i 

0  1  0  ±  sin  q^ 

0  0  1  0 

0  0  0  1 


M 


m,c, 


M 


<7jmi 


cos^ci  “Sin^^^^i 
sin^^,  cos^^i 
0  0 

0  0 

10  0  o' 

0  10  En,^ 

0  0  1  -Xg^ 

0  0  0  1 


0  0 
0  0 
1  0 
0  1 


(1.19) 

(1.20) 
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^*1«1 


M„,= 


sinr^i  0  -cos;'^^  0 

0  10  0 

cosrm,  0 

0  0  0  1 


0 

0 


sin^i 

0 

0 

1  COS^^/'l 

0 

0 

0 

1 

0 

0 

0 

1 

matrix  M 

c,p 

cc 

The  upper  and  lower  signs  of  the  element  5^^  sin^,  in  matrix  correspond  to  generation  of  right-  and  left-hand 

pinions,  respectively.  Wherein  the  modified  roll  is  applied  in  the  process  of  generation,  the  rotation  angles  and 
of  the  pinion  and  cradle  are  related  as 


bj  2  b'l  2 
V  ^  j 


=  »»ic(v'c,  -CwX  -Dwl) 


(1.21) 


where  ^2  ^3  modified  roll  parameters  and  C  and  D  the  modified  roll  coefficients. 


The  derivative  of  function  V'lCV'c,)  ta^'^en  at  v'c,  =0  determines  the  so-called  ratio  of  roll  represented  in  equation 
(1.21)  by  hi  or  Wi^. 

Equation  of  Meshing.  The  pinion  tooth  surface  Ej  is  the  envelope  to  the  family  of  head-cutter  surfaces.  Surface  E, 
consists  of  two  surfaces  that  correspond  to  the  cutter  surfaces  E^"^  and  E^*^ ,  respectively.  The  modified  roll  is  applied 
in  the  process  of  generation.  The  equation  of  meshing  for  generation  of  surface  Ep  ^  is  represented  as 

where  n^"^  is  the  unit  normal  to  the  surface,  and  v^f'^  is  the  velocity  in  relative  motion.  These  vectors  are  represented 

OTj  ^1 

in  the  fixed  coordinate  system  as  follows: 

(1.23) 

(1.24) 

vr  =  K? 

Here 


m\ 

rt 

Lm,c, 

1 

Ik’ 

-{o„0„  xfc»2> 

r«.  =M„,,M,,pr<‘'>(.p,^p) 

> 

It 

'o  ' 

-En. 

^bJ 

*"ml  “ 

[cosr„, 

0 

sinrmj 

0  OT<,,(5/,,)r 

The  ratio  is  not  constant  since  modified  roll  is  applied  and  can  be  represented  as 

ry,  _d¥c\  __  1 _ .  = _ ! - 

dt  dt  /w,^(l-2C(/^,-3Dv'd) 

where  C  and  D  are  modified  roll  coefficients. 

Finally,  we  obtain  the  equations  for  pinion  tooth  surface  E^"^  as 

r/‘’)(5p,^p,r,)  =  M,p(v^,)r^>(5p,0p) 

//;>(5p,0p,^,)  =  O 


(1.25) 


(1.26) 

(1.27) 
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Similarly,  the  pinion  fillet  surface  may  be  represented  as 


The  geometric  model  of  pinion  tooth  surfaces  is  shown  in  Fig.  1.11. 


(1.28) 

(1.29) 


Figure  1.1 1 :  3D  geometric  model  of  pinion  tooth  surfaces. 


1.5  Procedures  of  Proposed  Design 

The  authors  have  developed  an  approach  for  the  design  of  formate  cut  spiral  bevel  gears  that  provides  a  stabilized 
bearing  contact  and  reduced  level  of  noise. 

The  bearing  contact  is  of  a  longitudinal  direction  that  permits  an  increase  in  the  contact  ratio,  avoids  edge  contact,  and 
reduces  the  shift  of  the  bearing  contact  caused  by  misalignment. 

The  reduction  of  noise  is  achieved  by  application  of  a  predesigned  parabolic  function  and  limitation  of  maximal 
transmission  errors  caused  by  misalignment. 

The  approach  is  based  on  four  procedures  that  require  simultaneous  application  of  the  computational  algorithms 
developed  for  local  synthesis  and  simulation  of  meshing  and  contact  of  the  gear  drive. 

The  purpose  of  the  local  synthesis  is  to  determine  the  pinion  machine-tool  settings  for  the  following  conditions  of 
meshing  and  contact  of  pinion-gear  tooth  surfaces  Z,  and  Z2  ^ 

(i)  The  gear  machine-tool  settings  are  considered  as  given  (they  may  be  adapted,  for  instance,  from  a  Gleason 
summary  for  gear  generation). 

(ii)  The  pinion-gear  tooth  surfaces  Z,  and  Zj  are  in  tangency  at  a  chosen  mean  contact  point  M  (Fig.  1.12). 

(iii)  Parameters  2a ,  and  w'12  are  preassigned  and  determine  the  conditions  of  meshing  and  contact  of  Z| 

and  Z2  at  point  M  and  in  the  neighborhood  at  M  .  Here,  2a  is  the  major  axis  of  the  instantaneous  contaet 
ellipse,  TJ2  determines  the  tangent  to  the  path  of  contact  on  the  gear  tooth  surface  at  M ,  and  m\2  is  the 
derivative  at  point  M  where  m^2^,4>\)  is  the  gear  ratio. 

Note:  It  will  be  shown  in  Section  1.6  that  the  location  of  point  M2  on  the  gear  tooth  surface,  as  the  candidate  for  the 
mean  contact  point  of  Z|  and  Z2 ,  may  be  chosen  ahead.  Then,  the  procedure  of  the  local  synthesis  permits  determining 
such  a  point  M^  on  pinion  tooth  surface  Z,  that  will  be  in  tangency  with  Z2  at  the  chosen  point  M2  of  Zj. 
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Figure  1.12:  Illustration  of  parameters  %  and  a  applied  for  local  synthesis. 

The  simulation  of  meshing  and  contact  of  pinion-gear  tooth  surfaces  is  accomplished  by  application  of  the  developed 
TCA  (Tooth  Contact  Analysis)  computer  program.  The  inputs  for  TCA  are  the  equations  of  pinion-gear  tooth  surfaces, 
parameters  of  motion  and  assembly.  The  outputs  are  the  bearing  contact  and  the  function  of  transmission  errors. 

The  design  is  based  on  simultaneous  application  of  the  computational  algorithms  developed  for  local  synthesis  and 
TCA,  it  is  an  iterative  process,  and  requires  application  of  the  follovidng  four  procedures: 

Procedure  1:  The  path  of  contact  on  gear  tooth  surface  is  a  spatial  curve  L .  Projection  of  L  on  tangent  plane  T  (Fig. 
1.13)  is  designated  by  Xj- .  The  purpose  of  procedure  1  is  to  obtain  Lj  that  is  a  straight  line  directed  longitudinally. 

Figure  1.13  shows  projections  and  lip  that  deviate  from  the  desired  line  Lf . 


Figure  1.13:  Projections  of  various  paths  of  contact  Lp  on  tangent  plane  T . 
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Figure  1.14  shows  the  flow  chart  that  describes  the  steps  of  procedure  1.  The  procedure  of  computation  is  an  iterative 
process  and  requires  variation  of  parameter  w'j2 .  The  steps  of  the  procedure  are  as  follows: 


Figure  1.14:  Flow  chart  for  procedure  1:  design  of  bearing  contact. 

Step  1:  Input  the  gear  machine-tool  settings  into  the  program  of  local  synthesis. 

Step  2:  Input  parameters  2a,  7/2  » ^'12  ^he  program  of  local  synthesis. 

Step  3:  Apply  the  program  of  local  synthesis. 

Step  4:  Determine  from  the  computational  algorithm  of  local  synthesis  the  pinion  machine-tool  settings. 

Step  5:  Using  the  machine-tool  settings  obtained,  compute  the  pinion  and  gear  tooth  surfaces  Ej  and  E2  • 

Step  6:  Apply  TCA  computer  program  and  obtain  projection  Lj  of  the  path  of  contact  L  on  tangent  plane  T  (Fig. 

1.13). 

Step  7:  Apply  the  iterative  process  and  obtain  the  projection  of  the  path  of  contact  as  a  straight  line  by  variation  of 

Step  8:  Obtain  as  the  output  the  final  machine-tool  settings  for  the  pinion. 

The  iterative  process  is  based  on  the  following  considerations: 

(i)  The  TCA  computer  program  determines  the  projection  Lp  of  the  path  of  contact  numerically. 

(ii)  Consider  Lj  in  coordinate  system  Sj  (Fig.  1.13)  with  the  origin  coinciding  with  M  and  axes  (x^,y,)  are 
located  in  plane  T .  Axis  is  the  tangent  to  Lp  at  point  M . 

(iii)  Represent  projection  Ly  as  a  parabolic  regression  equation  by  using  a  regression  subroutine  [14] 

=  Po(m\2 )+ j5,(to’,2  +  )x^  (1.30) 
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where  /m',2  =  is  the  input  parameter  that  is  varied  in  the  process  of  iteration. 

(iv)  The  goal  is  to  represent  as  a  straight  line  that  might  be  obtained  by  such  a  magnitude  of  /n'jj  whereas 
coefficient  ^2  ('”'12  )  =  0  iterations  for  determination  of  >82  ("*'12  )  =  0  are  executed  by 

applying  the  secant  method  [13]  and  illustrated  by  Fig.  1.15.  Equation  (1.30)  with  ^2(^'i2)  ~  ®  provides  the 
projection  of  path  of  contact  as  a  straight  line. 


Figiue  1.15:  Illustration  of  computations  for  determination  of  )  =  0 

Procedure  2:  The  previously  discussed  procedure  1  provides  the  desired  longitudinal  path  of  contact  and  bearing 
contact.  However,  the  shape  of  the  obtained  function  of  transmission  errors  and  the  magnimde  of  maximal 
transmission  errors  do  not  satisfy  the  requirements  of  low-noise  gear  drive.  The  goal  of  procedure  2  is  to  obtain  a 
parabolic  function  of  negative  transmission  errors  and  of  limited  value  of  maximal  errors. 

The  sequence  of  steps  applied  for  procedure  2  is  represented  by  the  flow  chart  shown  in  Fig.  1.16. 

At  the  start,  we  use  the  pinion  and  gear  machine  tool  settings  obtained  in  procedure  1,  but  consider  that  the  pinion  is 
generated  by  application  of  modified  roll.  Then,  we  compute  the  gear  and  pinion  tooth  surfaces  S,  and  S2 »  use 
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them  for  the  TCA  computer  program.  The  intermediate  output  of  TCA  is  a  function  of  transmission  errors.  We 
compare  it  with  the  sought-for  parabolic  function  of  transmission  errors  and  determine  the  necessary  corrections  of 
modified  roll.  Repeating  the  procedure  discussed  above,  we  can  obtain  finally  the  sought-for  parabolic  function  of 
transmission  errors. 


(  start  ) 


_ ^ 

t 

1 

Gear  machine 
tool  settings 

Pinion  machine 
tool  settings 

Modified  roll 

r-  . :=zr 


Gear  tooth  surface 

Pinion  tooth  surface 

T 

TCA:  Simulation  of 
meshing  and  contact 


Output  predesigned 
parabolic  function  of 
transmission  errors 


C  End  3 


Figure  1,16:  Flow  chart  of  procedure  2, 

Analytically,  the  algorithm  of  procedure  2  is  represented  as  follows. 

Step  1:  Function  of  transmission  errors  A^(^)  obtained  from  TCA  is  numerically  represented  as  a  polynomial 
ftinction  up  to  third  member  to  be  included 

)  =  flo  +  (1.31) 

The  shape  and  magnitude  of  maximal  errors  of  function  do  not  satisfy  the  requirements  of  a  low-noise  gear 

drive.  Function  (1.31)  has  to  be  transformed  into  a  parabolic  function  with  limited  magnitude  of  transmission  errors  by 
application  of  procedure  2. 

Step  2:  We  apply  for  transformation  of  function  A^^^*^(^)  the  modified  roll  for  generation  of  the  pinion  that  is 
represented  as 

{¥c\ )  =  -  Vd  ( 1 .32) 

Here  y/^  and  are  the  angles  of  rotation  of  the  pinion  and  the  cradle  performed  during  the  process  for  generation, 
ftixc  is  the  first  derivative  of  y/\{y/c\)  taken  at  y/^^  obtained  by  the  procedure  of  local  synthesis  (see  Section 

1 .6).  Coefficients  62  and  have  to  be  determined  by  using  an  iterative  process  (see  below). 

Step  3:  Assigning  coefficients  bj  and  b^  of  modified  roll  and  using  pinion  machine-tool  settings  obtained  in 
procedure  I,  we  may  determine  the  pinion  tooth  surface  2j .  The  gear  tooth  surface  L2  Is  considered  as  given. 
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Step  4:  Applying  TCA  for  simulation  of  meshing  of  Ij  and  S2  ,  we  may  obtain  numerically  the  modified  function  of 
transmission  errors  and  then  represent  it  again  by  polynomial  function  (L3 1)  with  new  magnitudes  of  coefficients  , 
aj ,  a2  ,  and  . 

Procedure  2  with  steps  1  to  4  has  to  be  repeated  until  a  parabolic  function  of  transmission  errors  with  limited 

magnitude  of  maximal  transmission  errors  is  obtained. 

The  computation  is  an  iterative  process  based  on  application  of  secant  method  of  numerical  computations  [13].  The 
goal  is  to  transform  the  shape  of  obtained  function  )  of  transmission  errors  into  a  parabolic  function 

where 

( 

-j- 

This  goal  is  accomplished  by  variation  of  coefficients  f»2  ^  function  (1.32)  provided  by  modified  roll.  The 

variation  of  62  2nd  is  performed  independently  and  is  illustrated  by  drawings  of  Fig.  1.17. 


Figure  1.17:  Illustration  of  variation  of  coefficients  ^2  2*^^  ^  of  modified  roll. 
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Figure  1.17(a)  illustrates  variation  of  coefficient  of  modified  roll  to  obtain  coefficient  Oj  =  0  of  function  (1.31). 
Function  a2(^)  is  determined  as  the  output  of  TCA  by  variation  of  modified  roll. 

Figure  1.17(b)  illustrates  variation  of  coefficient  Z>2  of  modified  roll  for  obtaining  of  |a^(^)|^^^  =  A®  . 

Function  A^(^)  is  determined  as  the  output  of  TCA  by  variation  of  bj  • 

Procedure  3.  Procedures  1  and  2  are  performed  by  simultaneous  application  of  computational  algorithms  of  local 
synthesis  and  TCA,  but  the  errors  of  alignment  are  considered  at  these  procedures  as  equal  to  zero. 

The  purpose  of  procedure  3  is  to  adjust  the  gear  drive  to  the  existence  of  errors  of  alignment.  This  is  achieved  by 
adjusting  the  previously  designed  path  of  contact  to  the  assigned  or  expected  errors  of  alignment.  The  procedure  is 
based  on  the  flow  chart  shown  in  Fig.  1.18. 


Figure  1.18:  Flow  chart  of  procedure  3. 

The  adjustment  of  the  gear  drive  to  the  errors  of  alignment  is  performed  by  correction  of  parameter  72  (Fig-  1-12)  that 
determines  the  orientation  of  the  path  of  contact  on  the  gear  tooth  surface. 

The  prescribed  procedure  yields  that  in  some  cases  of  design  it  becomes  necessary  to  deviate  the  path  of  contact  from 
the  longitudinal  direction  to  reduce  the  shift  of  bearing  contact  caused  by  errors  of  alignment.  The  procedure  has  to  be 
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applied  for  separate  simulation  of  each  error  of  alignment  for  obtaining  the  required  correction  of  t]2  •  The  permissible 
tolerance  of  each  alignment  error  is  obtained  by  procedure  3. 

Procedure  4.  The  purposes  of  procedure  4  are  investigation  of  formation  of  bearing  contact  and  determination  of 
contact  and  bending  stresses  for  more  than  one  cycle  of  meshing.  The  goals  are  obtained  by  application  of  the  finite 
element  method  by  the  general  purpose  finite  element  computer  program  ABAQUS. 

Investigation  of  formation  of  bearing  contact  enables  to  discover  hidden  areas  of  severe  contact  due  to  the  elastic 
deformation  of  the  gear  teeth.  Such  a  contact  may  be  discovered  if  finite  element  models  of  several  pairs  of  teeth  are 
developed  and  analyzed  in  contact  positions  corresponding  to  more  than  one  cycle  of  meshing.  Such  contact  positions 
are  obtained  by  application  of  TC A  computer  program. 

Hidden  areas  of  severe  contact  are  accompanied  with  a  substantial  increase  of  contact  stresses.  Those  areas  of  severe 
contact  might  be  avoided  by  increasing  the  mismatch  of  generating  surfaces.  We  could  obtain  such  goal  by  application 
of  combination  of  straight-line  profile  blades  with  parabolic  profile  blades  of  the  pair  of  head-cutters  that  generate  the 
pinion  and  the  gear,  respectively. 

1.6  Local  Synthesis  and  Determination  of  Pinion  Machine  Tool  Settings 

As  a  reminder,  the  local  synthesis  is  accomplished  simultaneously  with  the  TCA,  and  the  computation  is  an  iterative 
process.  The  computer  program  developed  permits  the  determination  of  the  pinion  machine-tool  settings  at  each 
iteration. 

The  input  parameters  for  the  local  synthesis  are  the  gear  machine-tool  settings  and  parameters  7/2 ,  a ,  and  7w'i2  (see 
flow  chart,  Fig.  1.14). 

The  procedure  of  the  local  synthesis  is  as  follows: 

Stage  1:  Determination  of  the  mean  contact  point  on  the  gear  tooth  surface 

Ste/J  7;  Mean  point  A  on  surface  Z2  is  chosen  by  designation  of  parameters  and  (Fig.  1.19),  where  A  is  the 
candidate  for  the  mean  contact  point  M  of  surfaces  S2  •  Then  we  obtain  the  following  system  of  two 

equations  in  two  unknowns 

Z2(«;,<)  =  L.  (1-35) 

where  X2,  T2.and  Z2  are  the  projections  of  position  vector  r2(M*,^g). 
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Step  2:  Equations  (1.35)  considered  simultaneously  allow  us  to  determine  surface  parameters  for  point  A. 

Vector  fiinctions  and  ^g{0^)  determine  the  position  vector  and  surface  unit  normal  for  a  current  point  of 

gear  tooth  surface  .  Taking  in  these  vector  functions  =  w*  and  Bg^6\.  we  can  determine  the  position  vector 
of  point  A  and  the  surface  unit  normal  at  A . 


Step  3:  Parameters  Ug  and  0^  and  the  unit  vectors  and  of  principal  directions  on  surface  are  considered  as 
known.  Here: 


ar<^> 

-  cosOg  +  2a  ^Ug  sinotg 
(±  sin  OTg  ±  2a^Ug  cos  )  sin 

ds^ 

8 

(IsinOg  ±2a^«gCOsag)cos^g 

e 


u 


0 

8 

= 

cos^. 

de^ 

8 

g 

-sin^^ 

(1.36) 


(1.37) 


The  gear  generating  surface  is  a  surface  of  revolution,  and  the  principal  curvatures  k.  and  of  are  determined 

o  ^  S 

by  the  following  equations: 


=■ 


2o. 


0  +  4ayj 


*„  = 


-  cosag  4*  la^Ug  sina^ 


±s 


gosinag  ± Ug sinor^  ± a^Ug cosUg 


yi+4^ 


The  upper  and  lower  signs  correspond  to  the  concave  and  convex  sides  of  gear  tooth,  respectively. 


(1.38) 


The  unit  vectors  and  of  principal  directions  on  surface  Z2  are  obtained  by  using  matrix  transformation  from 
*^2  •  principal  curvatures  and  of  Z2  are  the  same  as  the  principal  curvatures  kg  and  , 

respectively,  of  the  gear  head-cutter. 


Stage  2:  Tangency  of  Surfaces  Z2(^g)  and  I,  at  Mean  Contact  Point  A/. 


Step  1:  The  derivations  accomplished  at  stage  1  permit  the  position  vector  and  the  surface  unit  normal  of 
point  A  of  tangency  of  surfaces  Z2(I^)  to  be  determined.  The  goal  now  is  to  determine  such  a  point  M  in  the  fixed 
coordinate  system  Si  (Fig.  1.20)  where  two  surfaces  X2(^g)  and  Z|  are  in  tangency  with  each  other. 


Surfaces  Z2  and  2^  are  the  same  ones,  because  there  is  no  relative  motion  between  the  cradle  and  workpiece 
(between  coordinate  system  ^2  to  Sg  )  whereas  the  gear  is  generated.  Using  coordinate  transformation  from  ^2  to  5/ 

(Fig.  1.20),  we  may  obtain  and  The  new  position  of  point  A  in  Si  will  become  point  M  of  tangency  of 
I2  and  Zj ,  if  the  following  equation  of  meshing  between  Z2  and  Zj  is  observed 

(1.39) 

Here,  =  and  s  is  the  relative  velocity  at  point  A  determined  with  the  ideal  gear  ratio 


„(0)  _  ® 
/K21'  =  — 
0) 


,(2) 


(') 


(1.40) 


The  solution  of  equation  (1 .39)  for  provides  the  value  of  the  turning  angle  .  It  is  evident  that  two  surfaces  Z2 

and  I,  are  now  in  tangency  with  each  other  at  point  M  . 


N  AS  A/CR— 2003 -212336 


24 


Figure  1.20:  Coordinate  systems  S2,  Si,  and  S',  applied  for  local  synthesis. 

Step  2:  Consider  at  point  M  that  the  principal  curvatures  and  kg  of  surface  Z2 ,  and  the  unit  vectors  and  of 
principal  directions  on  S2  know.  The  unit  vectors  and  are  represented  in  S, .  The  goal  now  is  to 

determine  at  point  M  the  principal  curvatures  kf  and  k,,  of  surface  S, ,  and  the  unit  vectors  and  e*  of  principal 
directions  on  S, .  This  goal  can  be  achieved  by  application  of  the  procedure  described  in  Section  1.7.  It  is  shown  in 
Section  1.7  that  the  determination  of  kf,  k,,,  ef  and  e*  becomes  possible  if  parameters  to’,2  ,  Vi  (or  ).  and  a!  S 
are  assumed  to  be  known  or  are  used  as  input  data. 

Stage  3:  Tangency  of  Surfaces  22 (^g)  >  '^p  Contact  Point  M . 

Basic  Equations:  Tangency  of  2)2  (2g)  and  S,  at  mean  contact  point  M  has  been  already  provided  at  the  previous 
stages.  The  position  vector  of  point  A/  and  the  surface  unit  normal  n^  ^  at  point  M  were  determined  in 
coordinate  system  S, .  Let  imagine  now  that  coordinate  system  5,  (it  coincides  with  S, )  and  surface  E,  are  installed 
in  coordinate  system  (Fig.  1.8).  Let  angle  (it  is  the  initial  value  of  )  be  the  installment  angle  of  the  pinion. 
Using  coordinate  transformation  from  5,  to  ,  we  may  determine  in  S„^  position  vector  of  point  M  and  the 
surface  unit  normal  .  The  conditions  of  improved  meshing  and  contact  of  pinion  and  gear  tooth  surfeces  E,  and 

/W] 

E2  are  considered  in  Section  1 .7  and  then  the  relationships  between  the  principal  curvatures  and  directions  of  surfaces 
for  such  conditions  of  meshing  and  contact  are  determined  (see  Eqs.  (1.81)).  The  point  of  tangency  of  siufaces  E,  and 
Ep  is  designated  in  Section  1.7  as  point  B  .  The  pinion  generating  surface  E^  is  installed  in  taking  the  cradle 
angle  v'c,  to  zero.  The  position  vector  of  point  B  of  surface  E^  and  the  surface  unit  normal  at  B  are 

represented  in  5^^  as  and  .  The  tangency  of  E,  and  E^  at  the  mean  contact  point  Af  is  satisfied,  if  the 
following  vector  equations  are  observed 


(1.41) 

Wj  mi 

S  : 

II 

(1.42) 

mi  mi 

.(A/)  (Ip)^O 

mi  mi 

(1.43) 
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where  equation  (1.43)  is  the  equation  of  meshing.  Observation  of  equations  (1.41)  to  (1.43)  means  that  all  of  the  three 
surfaces  (S2(Ig) ,  S|  and  )  are  in  contact  at  point  M . 


Using  equations  (1.41)  to  (1.43)  and  (1.81),  it  becomes  possible  to  obtain  the  settings  of  the  pinion  and  the  head-cutter 
that  guarantee  the  improved  conditions  of  meshing  and  contact  at  point  M .  The  machine  tool  settings  to  be 
determined  are  as  follows: 

(i)  ®id  Settings  Xg^  and  Xjj^  are  related  by  the  following  equation  (Fig. 

1..21) 

=  -(^£),  - 1  (1 .44) 

where 


M  °  ‘"'''1-  k  .  F./2)cos„ 


+  aJsmy^\ 


Equation  (1.44)  is  determined  for  the  case  of  design  of  a  spiral  bevel  gear  with  different  apexes  of  the  pitch 
and  root  cones. 


Figure  1.21:  For  derivation  of  pinion  machine  tool  settings. 

(ii)  Design  parameter  of  the  head-cutter  (Fig.  1.10). 

(iii)  Parameters  and  that  determine  the  installment  of  the  head-cutter  on  the  cradle  (Fig.  1.9). 


(iv)  Parameter  that  is  the  initial  angle  for  installment  of  coordinate  system  Sj  with  respect  to  (Fig. 
1.8),  and  parameter  0^  of  the  head-cutter  surface  . 

Computation  of  Pinion  Machine-Tool  Settings.  The  procedure  for  computation  is  as  follows; 


Step  I:  Determine  the  values  of  Op  and  that  are  the  sought-for  two  unknowns  in  Equations  (1.41)-(1.43). 
Equation  (1 .41)  is  used  for  determination  of  Op  and  ,  taking  into  account  that 


and 


cosOp  cos  Op 

cosUpSmOp 

Tsina^ 


(1.45) 


nW  =  L„.,L,,,n{ 

where  n{*^*  snj*^^  since  5,  coincides  with  S/  (Fig.  1.20).  Here, 


(1.46) 
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(1.47) 


0 

COS/m, 

0 

1 

0 

0 

sin/m. 

COS^I®^ 

-sin^^l^^  ( 

cos^[^^  ( 

Lm=, 

0  0  1 

Equations  (1.41)  and  (1.45)-(1.48)  yield  the  following  expressions  for  determination  of  0p  and 

Kfelsiny  sina 

cos  00  = - ‘ - 

cosy„j  cosap 

-%cos6'„cosap  +  %siny„,  sin^pCosap  +  n,^nfeCOS7'„j 

sin</f^  = — - - ; - f-2 - 51; - 

---.k+4) 


(0) 


.(0) 


sin^^  cosa^  -  % 

^ly 


(1.48) 

(1.49) 

(1.50) 

(1.51) 


where  is  the  given  value  of  the  profile  angle  of  the  head-cutter,  and  ,  n^y ,  and  rtj^  are  the  three  components  of 
vector  .  The  great  advantage  of  the  approach  developed  and  discussed  above  is  that  the  requirement  of  the 
coincidence  of  the  normals  does  not  require  a  non-standard  profile  angle  or  the  tilt  of  the  head-cutter  with  respect 
to  the  cradle. 

Using  0p ,  it  becomes  possible  to  determine  the  unit  vectors  and  of  principal  directions  on  surface  Zp  at  point 

M. 


5/^/?  2;  Determination  of  pinion  machine-tool  settings  m^p  and  the  design  parameter  Rp  of  the  head- 

cutter  (five  unknowns  are  sought-for). 

As  a  reminder,  X^^  and  are  related  by  equation  (1 .44).  The  determination  of  the  machine-tool  settings  mentioned 

above  is  based  on  application  of  the  system  of  equations  (1.81)  and  equation  (1.43)  that  represent  a  system  of  four  non¬ 
linear  equations  with  four  unknowns:  Xj)^ ,  ,  mip  and  Rp  . 


The  design  parameters  mentioned  above  provide  as  well  improved  conditions  of  meshing  and  contact  at  the  mean 
contact  point  M . 

Step  3:  Determination  of  machine-tool  settings  and  qi  and  the  pinion  surface  parameter  Sp  (three  xmknowns  are 
sought-for). 

Determination  of  the  mentioned  three  parameters  is  based  on  application  of  equation  (1.42),  considering  that 
generating  surface  is  a  cone.  The  final  equations  for  the  pinion  are  as  follows: 


cos^i  +  (Rp  +  Sp  smap)cos0p  = 

(1.52) 

sin{±gi)+(Rp  +SpSmap)sm0p 

(1.53) 

ry(M) 

-SpCOsap  =  Zl^> 

(1.54) 

The  upper  and  lower  signs  in  equation  (1.53)  correspond  to  the  design  of  right-hand  and  left-hand  pinions, 
respectively. 
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The  stages  of  computation  discussed  above  may  be  summarized  as  follows: 

(i)  It  is  necessary  to  determine  ten  unknowns:  six  machine-tool  settings  two 

surface  parameters  [Op.Sp),  one  cutter  parameter  and  one  position  parameter  which  defines  the 
pinion  initial  turn  angle, 

(ii)  The  equation  system  for  determination  of  the  unknovvms  is  formed  as  follows: 


In  addition,  we  use  three  curvature  equations 

tan2cT^’^^  = 


^33 


rtj] 

(1.55) 

JM)  ^  JB) 

(1.56) 

(1.57) 

(1.58) 

3^/23 

^23  V33 

‘2^/i3^23 

3sin2o’^^^^ 

(1.59) 

^33 

Equation  (1.55)  is  equivalent  to  two  independent  scalar  equations,  equations  (1.56)  are  equivalent  to  three  scalar 
equations,  equations  (1.57),  (1.58),  and  (1.59)  represent  five  scalar  equations.  Thus,  the  system  of  equations  provides 
indeed  ten  scalar  equations  for  determination  of  ten  unknowns.  The  solution  for  the  unknowns  requires:  (i)  solution  of 
a  subsystem  of  four  nonlinear  equations  (see  step  2),  and  (ii)  solution  of  six  remaining  equations  represented  in  echelon 
form  (each  of  the  six  equations  contains  one  unknown  to  be  determined). 


1.7  Relationships  Between  Principal  Curvatures  and  Directions  of  Mating 
Surfaces 


The  procedure  of  local  synthesis  requires  the  knowledge  of  principal  curvatures  and  directions  of  contacting  surfaces. 
In  the  case  of  design  discussed  above,  the  pinion  tooth  surface  is  represented  by  three  related  parameters  and  the 
determination  of  pinion  principal  curvatures  and  directions  is  a  complex  problem.  The  solution  to  the  problem  is 
simplified  by  presentation  of  pinion  principal  curvatures  and  directions  in  terms  of  principal  curvatures  and  directions 
of  the  head-cutter  and  parameters  of  motions  of  the  generating  process  (proposed  in  [2,  7]). 

Important  relations  between  principal  curvatures  and  directions  of  mating  surfaces  being  in  point  contact  (proposed  in 
[2, 7])  are  applied  for  local  synthesis  (see  Section  1 .6). 


Henceforth  two  types  of  instantaneous  contact  of  meshing  surfaces  will  be  considered:  (i)  along  a  line,  and  (ii)  at  a  point. 
Line  contact  is  provided  in  meshing  of  the  surface  being  generated  with  the  tool  surface.  Point  contact  is  provided  for 
meshing  of  the  generated  pinion  and  gear  tooth  surfaces.  The  determination  of  the  required  relationships  is  based  on  the 
approach  proposed  in  [4, 7].  The  basic  equations  in  the  approach  developed  are  as  follows  (see  the  nomenclature): 

(1.60) 

(1.61) 


=  h|,'l  +  X  n 


dt 


[n- V' 


(1.62) 


Equations  (1 .60)  and  (1 .61)  relate  the  velocities  of  the  contact  point  and  the  tip  of  the  unit  normal  in  their  motions  over 
the  contacting  surfaces.  Equation  (1.62)  represents  the  differentiated  equation  of  meshing.  Equations  (1.60)  to  (1.62) 
yield  a  skew-symmetric  system  of  three  linear  equations  in  two  unknowns  x,  and  of  the  following  structure 


(1  =  1-3) 


(1.63) 
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Here  ;ci  and  JCj  are  the  projections  of  the  velocity  of  the  contact  point  in  the  motion  over  one  of  the  surface  along  the 
principal  directions  of  the  mating  siudace.  In  case  of  line  contact  of  surfaces,  the  solution  for  the  unknowns  is  indefinite 
and  the  rank  of  the  system  matrix  of  the  linear  equations  is  one.  In  case  of  point  contact  of  surfaces,  the  solution  for  the 
unknowns  is  definite,  and  the  rank  of  the  system  matrix  is  two.  The  properties  mentioned  above  are  used  for  the 
derivation  of  the  sought  for  relationships  between  the  principal  curvatures  and  directions  of  the  meshing  surfaces. 


Meshing  of  Surfaces  Sj  ^i-  Surfaces  S,  and  Sj  are  in  point  contact  and  their  meshing  is  considered  in  fixed 
coordinate  system  S,  (Fig.  1.20).  Equations  (1.60)  to  (1.62)  yield  the  following  system  of  three  linear  equations  [4,  7] 

(/■  =  1  -  3)  (1  -64) 


where 

v®  =  v<‘>-e,  v<'>  =  v®-e^ 

aii=ks-kf  cos^ sin^ 

0x2  =  £^21  “  0.5(A:y^  -A:/,)sin2cr^*^^ 

«13  =  «31  =  +  [("  ^  ®  J 

a22  =  kg-kf  sin^  -  k/,  cos^ 

"23  =  «32  =  +  I"  ^ 

fl33  =  -  [{n  X 

-  n  •  X  v^^)-  X  v^^)]+  /k',2  (n  x  kj)-  r 


(1.65) 


(1.66) 


It  is  considered  as  known  in  the  following  considerations:  point  M  of  tangency  of  surfaces  Sj  and  S2 ,  the  common 
surface  unit  normal,  the  relative  velocity  ,  the  principal  curvatures  and  and  directions  and  on  23  at 
M ,  and  the  elastic  deformation  S  of  surfaces  at  M .  The  goal  is  to  determine  the  principal  curvatures  kf  and  k^  and 

the  angle  formed  by  the  unit  vectors  e j-  and  . 


The  velocity  (/  =  1,2)  of  the  contact  point  in  motion  over  surface  2^  has  a  definite  direction  and  therefore  the  rank 
of  the  system  matrix  (1.64)  is  one.  This  property  yields  the  following  relation: 

\^U  ^12  ^13 

=  F(kj^,kf,,k,,k^,a^^^\m\2)  =  0  (1*67) 


^21 

P31 


^22 

^32 


^23 

^33 


The  sought-for  solution  for  kf ,  k,,,  and  can  be  obtained  if  the  following  parameters  will  be  chosen:  the 

derivative  m'j2  ;  the  ratio  alS ,  where  a  is  the  major  axis  of  the  contact  ellipse;  the  direction  at  M  of  the  tangent  to 
the  contact  path  on  one  of  the  couple  of  contacting  surfaces  Zj  and  £2  •  The  relation  between  the  directions  at  M  of 
the  tangents  to  the  contact  paths  on  both  surfaces  is  represented  by  the  equation  [4,  7] 


Xsarjx  =- 


'(‘^33 


,(12) 


+  «13V^ 


)tan;72 


O33  +  022 


(1.68) 


Choosing  772  at  point  M ,  we  can  determine  r]\ . 
Procedure  of  Determination  of  kf,  k,,  and 
5/e/»  J:  Determine  %  choosing  r\2. 
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Step  2  : 


(1.69) 


(1.70) 

(1.71) 

(1.72) 

(1.73) 

(1.74) 

(1.75) 


Step  8  : 

where  =k^+k^ . 

(1.76) 

Step  9  : 

*/  =  te^  +  gi)/2 

(1.77) 

Step  JO: 

(1.78) 

The  procedure  represented  above  can  be  used  to  obtain  the  sought  for  principal  curvatures  and  at  point  M  of 
tangency  of  surfaces  1,2  and  S,  and  the  principal  directions  on  Z,  at  A/ . 

Meshing  of  Surfaces  I,  and  :  The  tool  surface  Z^  generates  the  pinion  tooth  surface  Z, .  Surfaces  Z^  and  Z, 
are  in  line  contact  and  point  B  is  the  chosen  point  of  the  instantaneous  line  of  contact.  The  meshing  of  surfaces  is 
considered  in  .  At  point  B  ,  the  following  parameters  are  assumed  as  known:  the  curvatures  ky  and  k^  of  surface 

Z| ;  the  unit  vectors  Oy  and  of  principal  directions  on  Z| ;  the  surfaces  unit  normals;  the  relative  velocity  . 

The  goal  is  to  determine  the  principal  curvatures  kp  and  k,  of  surface  Zp  ,  and  the  angle  that  is  formed  by  the 
unit  vectors  Cy  and  Cp. 

Equation  (1.60)  to  (1.62)  yield  a  system  of  three  linear  equations 

=  ('  =  >-3)  (1.79) 
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where 


V^P'>  =  -Cy  , 


(1,80) 


The  direction  of  is  indefinite  since  surfaces  and  Zj  are  in  line  contact.  Therefore,  the  rank  of  system  matrix 
of  equations  is  equal  to  one.  Using  this  property,  the  following  equations  are  obtained: 

~2^13<^23 _ 


tan2<T^^^^  =- 


^23  ■“  <^13  ~  (^/  “  ^/|)^33 


k.  -k 

f  J 


^  rf33sin2cr^^^^ 


(1.81) 


Here 


where 


k,+kp  =  kf  +  k^+^^^^^ 

"33 

</i3  =  +  [(n  X  e  yj 

<^23  =  +  [(n  X  ej 

-  n  •  [(w^'^  X  X  v^'^)] 


e* 


(1.82) 


(1.83) 


1.8  Simulation  of  Meshing  and  Contact  Accomplished  in  Procedure  1  (Section  1.5) 

The  main  goal  of  simulation  of  meshing  and  contact  in  procedme  1  (see  flow  chart  of  Fig.  1.14)  is  the  determination  of 
bearing  contact  that  corresponds  to  pinion  machine-tool  settings  obtained  in  Section  1 .6.  Recall  that  the  computational 
procedures  represented  in  Sections  1.5  and  1.6  must  be  applied  simultaneously,  to  represent  an  iterative  process  for 
determination  of  longitudinal  path  of  contact.  The  developed  TCA  (Tooth  Contact  Analysis)  computer  program  (based 
on  simulation  of  meshing  and  contact)  allows  us  to  obtain  at  each  stage  of  iteration  the  path  of  contact.  In  addition, 
TCA  enables  us  to  obtain  the  function  of  transmission  errors  which  shape  has  to  be  corrected  by  application  of  pinion 
modified  roll  (see  procedure  2  of  Section  1 .5). 

Applied  Coordinate  Systems.  The  meshing  of  gear  tooth  surfaces  is  considered  in  the  fixed  coordinate  system  5* 
that  is  rigidly  connected  to  the  housing  (Fig.  1.22).  Movable  coordinate  system  Sj  and  S2  are  rigidly  connected  to  the 
pinion  and  the  gear,  respectively.  Auxiliary  coordinate  systems  and  are  used  for  representation  of  the  rotation 
of  the  pinion  (with  respect  to  )  and  the  gear  (with  respect  to  ).  The  errors  of  alignment  are  simulated  by 
respective  installment  of  and  5*^  with  respect  to  5^ . 


The  errors  of  installment  are:  A4,  -the  axial  displacement  of  the  pinion;  Ay  -the  change  of  the  shaft  angle  y ;  AE  -the 
shortest  distance  between  the  axes  of  the  pinion  and  the  gear  when  the  pinion-gear  axes  are  crossed  instead  of 
intersected;  A/42-the  axial  displacement  of  the  gear.  In  the  case  of  an  aligned  gear  drive,  we  consider  that  A4, ,  Ay, 

AE ,  and  AA2  are  all  equal  to  zero. 

Simulation  Algorithm.  During  the  process  of  meshing,  the  pinion  and  gear  tooth  surfaces  must  be  in  continuous 
tangency  and  this  condition  can  be  provided  if  at  any  instant  their  position  vectors  coincide  and  the  surface  normals  are 
collinear.  Instead  of  collinearity  of  surface  normals,  equality  of  surface  unit  normals  may  be  required. 
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Figure  1 .22;  Coordinate  systems  applied  for  simulation  of  meshing  of  pinion  and  gear  tooth  surfaces. 


Pinion  and  gear  tooth  surfaces  are  represented  in  coordinate  system  S/,  by  the  following  equations 

rh\^p^^p,WiA)  =  M**,  ,  (^  )r,  (Sp  ,0p,W\) 

Here 


cos^  sin^  0  0 


M 


i,i  - 


M 


AA,  - 


-sin^  cos^  0  0 
0  0  10 
0  0  0  1 
1  0  0  0  ■ 

0  10  0 

0  0  1  A/4i 

0  0  0  1 


M 


ba  - 


cos^ 

sin^ 

0 

0 


-sin^ 

cos^ 

0 

0 


0  0 
0  0 
1  0 
0  1 


cos(y  +  Ay)  0 
0  1 

sin(y  +  Ay)  0 
0  0 


-  sin(y  +  Ay) 
0 

cos(y  +  Ay) 
0 


-  A/42  sin(y  +  Ay) 
AE 

A^2  cos(y  +  Ay) 

1 


(1.84) 

(1.85) 
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Unit  normals  to  surfaces  of  pinion  and  gear  are  represented  in  by  the  following  equations,  respectively 


Here; 


Lai  - 


cos^ 

-sin^ 

0 


sin^i  0 
cos^  0 
0  1 


cos^  -sin^  0 
sin^  cos^2  ® 

0  0  1. 
cos(y  +  Ay)  0  -  sin(y  +  Ay) 
0  1  0 
sin(y  +  Ay)  0  cos(y  +  Ay) 


(1.86) 

(1.87) 


Conditions  of  continuous  tangency  of  pinion  and  gear  tooth  surfaces  are  represented  by  the  following  equations 


rf\sp,ep,\i/i,(l\)-rl^\ug,eg,^)  =  0 

(1.88) 

(1.89) 

fl;\sp,ep,ii7,)^Q 

(1.90) 

Surfaces  Sj  and  Ej  are  represented  in  Sf,  by  four  and  three  related  parameters,  respectively.  Equation  (1.90)  is  the 
equation  of  meshing  of  the  pinion  and  generating  head-cutter  and  relates  the  surface  parameters  of  the  head-cutter  with 
the  generalized  parameter  of  motion  of  the  process  of  generation  of  the  pinion.  Vector  equations  (1 .88),  (1 .89),  and 

(1.90)  yield  three,  two  independent  and  one  scalar  equations,  respectively.  The  whole  system  of  equations  (1.88)  to 
(1.90)  provides  six  equations  for  determination  of  seven  unknowns  represented  as 

fi{sp,0p,y/iA^Ug, dg,^)=0,  fi^C\  (i  =  l,-,6)  (1.91) 

One  of  the  unknowns,  say  ^ ,  is  considered  as  the  input  parameter  in  the  range 


71 


It  is  assumed  that  the  Jacobian  of  the  system  (1.88)  to  (1.90)  differs  from  zero  at  each  iteration. 


The  paths  of  contact  on  the  pinion  and  the  gear  tooth  surfaces  are  represented  by  the  following  functions: 


(1.92) 

(1.93) 

and 

(1.94) 

respectively. 

The  fiinction  of  transmission  errors  is  defined  as 

(1.95) 

The  bearing  contact  is  formed  as  a  set  of  instantaneous  contact  ellipses.  The  lengths  of  the  major  and  minor  axes  of 
contact  ellipse  and  their  orientation  are  determined  using  the  approach  proposed  in  [7]. 
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Using  vector  function  (1.94) ,  we  may  determine  projection  Lp  of  the  path  of  contact  on  the  gear  tooth  surface  in  the 
tangent  plane  at  the  mean  contact  point  M .  The  goal  is  to  obtain  as  a  straight  line  directed  longitudinally.  The 
process  of  computation  is  an  iterative  process  and  requires  simultaneous  application  of  local  synthesis  and  TCA.  The 
errors  of  alignment  in  this  procedure  are  not  taken  into  account.  The  influence  of  errors  of  alignment  on  transmission 
errors  requires  additional  application  of  TCA. 

It  may  be  discovered  in  some  cases  of  design,  that  the  gear  drive  is  too  sensitive  to  errors  of  alignment.  Then,  it 
becomes  necessary  to  deviate  Lp  from  the  longitudinal  directed  straight  line  by  correction  of  t]2  in  procedure  3 
(see  flow  chart  of  Fig.  1.1 8). 
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Chapter  2 

TCA  and  FEA  Results 


2.1  Introduction 

In  this  chapter,  the  results  obtained  from  computations  of  tooth  contact  analysis  (TCA)  and  finite  element  analysis 
(FEA)  are  represented.  Several  cases  of  the  design  of  formate-cut  spiral  bevel  gears  and  face-milled  generated  spiral 
bevel  gears  have  been  considered  for  their  comparison.  A  predesigned  fimction  of  transmission  errors  of  maximum 
level  of  8  arcsec  has  been  considered  for  all  cases  of  design. 

2.2  Design  Cases 

A  face-milled  generated  spiral  bevel  gear  drive  of  19  x  62  teeth  has  been  considered  as  the  baseline  design.  Our  goal  is 
to  find  an  optimized  design  of  this  gear  drive  using  the  formate  cut  method  providing  a  predesigned  fimction  of 
transmission  errors  of  low  level  and  a  stabilized  and  longitudinal  path  of  contact.  Using  the  formate-cut  method,  the 
productivity  of  manufacturing  of  the  gear  will  be  increased  substantially. 

The  following  cases  of  design  have  been  considered: 

Design  GTDl:  This  case  of  design  is  the  baseline  design  with  a  path  of  contact  directed  longitudinally.  Gear  and 
pinion  are  generated  by  head-cutters  with  straight  profiles. 

Design  GTDl:  In  this  case  of  design,  the  gear  is  generated  with  a  head-cutter  of  parabolic  profiles  whereas  the  pinion 
is  generated  by  head-cutters  with  straight  profiles. 

Design  GTD3:  In  this  case  of  design,  the  gear  is  generated  with  a  head-cutter  of  straight  profiles  whereas  the  pinion  is 
generated  by  head-cutters  with  parabolic  profiles. 

Design  FTDl:  In  this  case  of  design,  the  formate-cut  method  is  considered  for  the  gear.  The  head-cutter  for  the  gear 
has  parabolic  profiles.  The  pinion  is  generated  by  head-cutters  with  straight  profiles. 

Design  FTDl:  In  this  case  of  design,  the  formate-cut  method  is  considered  for  the  gear  as  well.  The  head-cutter  that 
cut  the  gear  has  parabolic  profiles,  but  the  parabola  coefficient  is  larger.  The  pinion  is  generated  by  head-cutters  with 
straight  profiles. 
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2.3  Results  of  Tooth  Contact  Analysis  (TCA) 

Only  results  of  tooth  contact  analysis  for  the  convex  side  of  the  gear  in  mesh  with  the  concave  side  of  the  pinion  will 
be  shown.  Similar  results  can  be  obtained  for  the  concave  side  of  the  gear  in  mesh  with  the  convex  side  of  the  pinion. 

Figure  2.1  shows  the  path  of  contact  and  the  function  of  transmission  errors  for  the  baseline  design.  Figures  2.2  and  2.3 
show  the  path  of  contact  and  the  function  of  transmission  errors  for  designs  GTD2  and  GTD3,  respectively.  Figures  2.4 
and  2.5  show  the  path  of  contact  and  the  function  of  transmission  errors  for  designs  FTDl  and  FTD2,  respectively, 
based  on  the  formate  cut  method. 


-10  000  -MDOO  -10  000  0  000  10  000  JO  000 

Pin.  Rotation  Angle  (dog! 


(b) 

Figure  2.1.  (a)  Path  of  contact  and  (b)  function  of  transmission  errors  for  case  of  design  GTDl, 
considered  as  well  as  the  baseline  design. 
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Trans. Eccot  (arc  sac)  Tcans. Error  (arc  sec) 


(a) 


Figure  2.5.  (a)  Path  of  contact  and  (b)  function  of  transmission  errors  for  FTD2  design. 
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2.4  Investigation  of  Limits  and  Tolerances  to  Errors  of  Alignment 

The  limits  and  tolerances  for  absorption  of  errors  of  alignment  have  been  investigated  for  the  previously  mentioned 
five  cases  of  design.  The  limits  of  the  values  of  errors  of  alignment  that  can  be  absorbed  by  the  gear  drive  are  those  that 
make  the  path  of  contact  for  one  cycle  of  meshing  to  have  one  contact  point  on  the  edge  of  the  gear  tooth  surface  (see 
Fig.  2.6).  Tables  2.1  to  2.5  show  the  limits  and  tolerances  to  errors  of  alignment  for  the  five  cases  of  design 
represented  above. 


Table  2.1.  Limits  and  to 

erances  to  errors  of  alignment  for  the  case  of  design  GTDl . 

Name  of  Parameters 

Pinion  concave-Gear  convex 

Min. 

Max. 

Tolerance 

Axial  Displacement  of  the  pinion, 
A/4j  (mm) 

-0.02  mm 

0.45  mm 

0.47  mm 

Axial  Displacement  of  the  gear, 

-0.23  mm 

0.55  mm 

0.78  mm 

Shortest  distance  between  axes  of 
pinion  and  gear,  AJ?  (mm) 

0.08  mm 

0.30  mm 

0.22  mm 

Change  of  crossing  angle,  Af 
tdeel 

-0.09  deg 

0.40  deg 

0.49  deg 

1 _ ^ 

Table  2.2.  Limits  and  tolerances  to  errors  of  alignment  for  the  case  of  design  GTD2. 

Name  of  Parameters 

Pinion  concave-Gear  convex  J 

Min. 

Max. 

Tolerance 

Axial  Displacement  of  the  pinion, 
AAj  (mm) 

-0.15  mm 

0.40  mm 

0.55  mm 

Axial  Displacement  of  the  gear, 
AA2  (mm) 

-0.15  mm 

0.90  mm 

1.05  mm 

Shortest  distance  between  axes  of 
pinion  and  gear,  AE  (mm) 

0.11  mm 

0.35  mm 

0.24  mm 

Change  of  crossing  angle,  A;' 
fdes) 

-0.13  deg 

0.36  deg 

0.49  deg 

Table  2.3.  Limits  and  tolerances  to  errors  of  alignment  for  the  case  of  design  GTD3. 

Name  of  Parameters 

Pinion  concave-Gear  convex 

Min. 

Max. 

Tolerance 

Axial  Displacement  of  the  pinion, 
A4j  (mm) 

-0.16  mm 

0.37  mm 

0.53  mm 

Axial  Displacement  of  the  gear, 
AA2  (mm) 

-0.11  mm 

0.68  mm 

0.79  mm 

Shortest  distance  between  axes  of 

pinion  and  gear,  AE  (mm) 

0.10  mm 

0.35  mm 

0.25  mm 

Change  of  crossing  angle.  A/ 

(deg) _ _ 

-0.16  deg 

0.33  deg 

0.49  deg 
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Table  2.4.  Limits  and  tolerances  to  errors  of  alignment  for  the  case  of  design  FTDl . 


Name  of  Parameters 

Pinion  concave-Gear  convex  | 

Min. 

Max. 

Tolerance 

Axial  Displacement  of  the  pinion, 
AAj  (mm) 

-0.05  mm 

0.06  mm 

0.11  mm 

Axial  Displacement  of  the  gear, 
A/42(mm) 

-0.05  mm 

0.10  mm 

0. 1 5  mm 

Shortest  distance  between  axes  of 
pinion  and  gear,  AE  (mm) 

-0.03  mm 

0.03  mm 

0.06  mm 

Change  of  crossing  angle,  Ay 
(deg) 

-0.09  deg 

0.07  deg 

0.1 6  deg 

Table  2.5.  Limits  and  tolerances  to  errors  of  alignment  for  the  case  of  design  FTD2. 


Name  of  Parameters 

Pinion  concave-Gear  convex  | 

Min. 

Max. 

Tolerance 

Axial  Displacement  of  the  pinion, 
AAj  (mm) 

-0.10  mm 

0.10  mm 

0.20  mm 

Axial  Displacement  of  the  gear, 
AA2  (mm) 

-0.08  mm 

0.10  mm 

0.18  mm 

Shortest  distance  between  axes  of 
pinion  and  gear,  AE  (mm) 

-0.04  mm 

0.05  mm 

0.09  mm 

Change  of  crossing  angle.  Ay 
(deg) 

-0.12  deg 

0.1 1  deg 

0.23  deg 

Figure  2.6.  Case  of  design  FTD2:  shift  of  the  bearing  contact  for  (a)  shortest  distance  between  axes  of  pinion  and  gear 
A£  =  -0.04  mm  and  (b)  shortest  distance  between  axes  of  pinion  and  gear  Af  =  0.05  mm. 

For  cases  of  design  GTDl,  GTD2,  and  GTD3,  the  following  errors  of  alignment  have  been  compensated  in  generation: 

Axial  displacement  of  the  pinion,  AA^  =  0.23  mm. 

Shortest  distance  between  axes  of  pinion  and  gear,  AE  =  0.16  mm. 

Change  of  crossing  angle,  Ay  =0.121  deg. 

No  errors  of  alignment  have  been  compensated  for  cases  of  design  FTDl  and  FTD2.  However,  this  capability  is 
available  in  the  developed  computer  program. 
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Figure  2.7  shows  a  graph  with  the  tolerances  to  errors  of  alignment  for  the  five  cases  of  design  considered.  The 
following  conclusions  can  be  drawn  fi'om  the  obtained  tolerances  to  errors  of  alignment. 

(i)  Formate-cut  spiral  bevel  gear  drives  (designs  FTDl  and  FTD2)  are  more  sensitive  to  all  the  different  types  of 
alignment  errors  than  generated  spiral  bevel  gear  drives  (designs  GTDl,  GTD2,  and  GTD3). 

(ii)  Spiral  bevel  gear  drives  are  very  sensitive  to  error  of  alignment  A£^ . 

(iii)  The  more  mismatch  between  the  profiles  of  the  head-cutters,  the  less  sensitivity  to  errors  of  alignment  (see 
tolerances  for  cases  FTDl  and  FTD2). 


TOLERANCES  TO  ERRORS  OF  ALIGNMENT 

1.2  n 


,  . . .in.  , - - - ^ ^ 

GTD  GTD  GTD  FTDl  FTD 


I H  Delta  A1  (mm)  H  Delta  A2  (mm)  ^  Delta  E  (mm)  B  Delta  gamma  (deg) _ 

Figure  2.7.  Tolerances  to  errors  of  alignment  for  five  cases  of  design  for  generated  and 
formate  cut  spiral  bevel  gear  drives. 

As  shown  in  Figure  2.7,  spiral  bevel  gear  drives  are  very  sensitive  to  the  change  of  the  shortest  distance  between  axes, 
AE .  A  method  to  compensate  the  shift  of  the  path  of  contact  caused  by  error  of  alignment  AE  is  proposed.  This 
method  is  based  on  compensation  of  error  of  aligiunent  AE  by  the  axial  displacement  of  the  pinion  AA^ .  Figure  2.8 
shows  an  example  of  compensation  of  error  of  alignment  error  of  alignment  AE  for  the  case  of  a  formate  cut  spiral 
bevel  gear  drive.  The  case  of  design  FTD2  has  been  considered  for  the  compensation  shown  in  Figure  2.8. 

Figure  2.8(a)  shows  the  path  of  contact  when  no  errors  of  alignment  occur.  Figure  2.8(b)  shows  the  path  of  contact  for 
an  error  of  alignment  AE  =0.04  nun.  Figure  2.8(c)  shows  the  path  of  contact  for  an  error  of  alignment  AE  =0.04  mm 
along  with  an  axial  displacement  of  the  pinion  A/4j=0.15  mm.  As  shown  in  Figure  2.8(c),  an  axial  displacement  of  the 
pinion  can  compensate  the  shift  of  the  path  of  contact  caused  by  the  shortest  distance  between  axes  AE .  The  same 
compensation  can  be  achieved  with  face-milled  generated  spiral  bevel  gears. 


NASA/CR— 2003-2 12336 


41 


Figure  2.8(d)  shows  the  function  of  transmission  errors  when  error  of  alignment  Af  =0.04  mm  is  compensated  with  an 
axial  displacement  of  the  pinion  AA^  =0.15  mm.  The  transmission  errors  are  of  parabolic  shape  with  a  maximum  level 
of  8.8  arcsec. 


Pin.  Rotation  Angle  (deg) 

Figure  2.8.  Case  of  design  FTD2:  (a)  path  of  contact  when  no  errors  of  alignment  occur,  (b)  path  of  contact  for  error  of 
alignment  AE  =0.04  mm,  (c)  path  of  contact  when  error  of  alignment  AE  =0.04  mm  is  compensated  by  A/4,  =0. 1 5 
mm,  (d)  function  of  transmission  errors  for  conditions  of  item  (c). 
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2.5  Enhanced  Approach  for  Finite  Element  Analysis  (FEA) 

Application  of  finite  element  analysis  permits: 

•  Determination  of  contact  and  bending  stresses  for  pinion  and  gear. 

•  Investigation  of  formation  of  bearing  contact  and  transfer  of  load  to  the  next  contacting  pair  of  teeth. 

•  Investigation  of  existence  of  hidden  areas  of  edge  contact  in  high-loaded  gear  drives. 

Finite  element  analysis  has  been  performed  by  application  of  the  general  purpose  computer  program  ABAQUS. 

An  enhanced  approach  for  performing  of  finite  element  analysis  that  has  the  following  advantages  has  been  developed: 

(a)  The  same  programming  language  (FORTRAN)  is  applied  for  synthesis,  analysis,  and  generation  of  finite 
element  models  of  the  gear  drives.  Graphical  interpretation  of  the  output  is  obtained  by  using  the  Visual 
Fortran  Quick  Win  run  time  library. 

(b)  The  generation  of  the  finite  element  mesh  required  for  finite  element  analysis  is  performed  automatically 
using  the  equations  of  the  surfaces  of  the  tooth  and  the  portion  of  the  rim.  Nodes  of  the  finite  element 
mesh  lying  on  the  tooth  surfaces  of  pinion  (gear)  are  guaranteed  to  be  points  of  the  real  tooth  surfaces  of 
the  pinion  (gear).  Loss  of  accuracy  due  to  the  development  of  solid  models  using  CAD  computer 
programs  is  avoided.  The  boundary  conditions  for  the  pinion  and  the  gear  are  set  automatically  as  well. 

(c)  Investigation  of  the  bearing  contact  formation  for  the  whole  cycle  of  meshing  permits  the  discovery  of 
hidden  areas  of  edge  contact. 

(d)  There  is  no  need  to  apply  CAD  computer  programs  for  the  development  of  finite  element  models 
required  for  application  of  general  purpose  FEA  computer  program  ABAQUS. 

The  development  of  the  solid  models  and  finite  element  meshes  using  CAD  computer  programs  is  time  expensive, 
requires  skilled  users  of  used  computer  programs  and  has  to  be  done  for  every  case  of  gear  geometry  and  position  of 
meshing  desired  to  be  investigated. 

The  developed  approach  is  fi-ee  of  all  these  disadvantages  and  is  summarized  as  follows: 

Step  1;  Using  the  equations  of  both  sides  of  tooth  surfaces  and  the  portions  of  the  corresponding  rim,  we  may 
represent  analytically  the  volume  of  the  designed  body.  Figure  2.9(a)  shows  the  designed  body  for  one-tooth  model  of 
the  pinion  of  a  spiral  bevel  gear  drive. 

Step  2:  Auxiliary  intermediate  surfaces  1  to  6  as  shown  in  Fig.  2.9(b)  can  be  determined.  Surfaces  1  to  6  enable  us  to 
divide  the  tooth  in  six  parts  and  control  the  discretization  of  these  tooth  subvolumes  into  finite  elements. 

Step  3:  Analytical  determination  of  node  coordinates  is  performed  taking  into  account  the  number  of  desired  elements 
in  longitudinal  and  profile  direction  (Fig.  2.9(c)).  We  emphasize  that  all  nodes  of  the  finite  element  mesh  are 
determined  analytically  and  those  lying  on  the  intermediate  siirfaces  of  the  tooth  are  indeed  points  belonging  to  the  real 
surface. 

Step  4:  Discretization  of  the  model  by  finite  elements  using  nodes  determined  in  previous  step  is  accomplished  as 
shown  in  Fig.  2.9(d). 
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Figure  2.9.  Illustrations  of  (a)  the  volume  of  designed  body,  (b)  auxiliary  intermediate  surfaces,  (c)  determination 
of  nodes  for  the  whole  volume,  and  (d)  discretization  of  the  volume  by  finite  elements. 


Step  5:  Setting  of  boundary  conditions  for  gear  and  pinion  is  performed  automatically.  Nodes  on  the  sides  and  bottom 
part  of  the  rim  portion  of  the  gear  are  considered  fixed  (Fig.  2.10(a)).  Nodes  on  the  two  sides  and  bottom  part  of  the 
rim  portion  of  the  pinion  build  a  rigid  surface  (Fig.  2.10(b)).  Rigid  surfaces  are  three-dimensional  geometric  structures 
that  cannot  be  deformed  but  can  perform  translation  or  rotation  as  rigid  bodies.  They  are  also  very  cost-effective  since 
the  variables  associated  with  a  rigid  surface  are  the  translations  and  rotations  of  a  single  node,  known  as  the  rigid  body 
reference  node  (Fig.  2.10(b)).  The  rigid  body  reference  node  is  located  on  the  pinion  axis  of  rotation  with  all  degrees  of 
freedom  but  the  rotation  around  the  axis  of  rotation  of  the  pinion  fixed  to  zero.  The  torque  is  applied  directly  to  the 
remaining  degree  of  freedom  of  the  rigid  body  reference  node  (Fig.  2.10(b)). 


Step  6:  Definition  of  contacting  surfaces  for  the  contact  algorithm  of  the  finite  element  computer  program  ABAQUS  is 
obtained  automatically  as  well  and  requires  definition  of  the  master  and  slave  surfaces.  Generally,  the  master  surface  is 
chosen  as  the  surface  of  the  stiffer  body  or  as  the  surface  with  the  coarser  mesh  if  the  two  surfaces  are  on  structures 
with  comparable  stiffness.  Finite  element  mesh  for  the  pinion  is  more  refined  than  the  one  for  the  gear  due  to  larger 
curvatures.  The  gear  and  pinion  tooth  surfaces  are  considered  as  the  master  and  slave  surfaces  respectively  for  the 
contact  algorithm. 
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Boundary  Conditions  for  the  Gear 


(b) 


Figure  2.10.  (a)  Boundary  conditions  for  the  gear  and  (b)  schematic  representation  of 
boundary  conditions  and  application  of  torque  for  the  pinion. 


2.6  Results  of  FEA 

Five  cases  of  design  of  face-milled  generated  spiral  bevel  gear  drives  and  formate  cut  spiral  bevel  gear  drives  have 
been  accomplished. 


The  discretization  of  pinion  and  gear  teeth  for  a  3-pair-of-teeth  finite  element  model  is  shown  in  Fig.  2.11. 
Incompatible  mode  elements  C3D8I  have  been  used  to  build  the  finite  element  mesh.  These  elements  are  first-order 
elements  that  are  enhanced  by  incompatible  modes  of  deformation  to  improve  their  bending  behavior.  The  total 
number  of  elements  is  33606  with  42076  nodes  for  each  geometric  model.  The  material  is  steel  with  the  properties  of 
Young's  Modulus  £’=2.068  10*  mN/mm^  and  the  Poisson's  ratio  0.29.  A  torque  of  549  N  m  has  been  applied  to  the 
pinion  of  the  gear  drives.  The  finite  element  mesh  of  the  whole  gear  drive  for  case  of  design  FTDl  is  shown  in 
Fig.  2.12. 
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Figure  2.11.  Finite  element  model  of  three  pairs  of  contacting  teeth  used  for  FEA. 


Figure  2.12.  Whole  finite  element  mesh  of  the  formate-cut  spiral  bevel  gear  drive  of  design  FTDl . 
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The  concave  side  of  the  pinion  and  the  convex  side  of  the  gear  are  considered  as  the  driving  and  driven  surfaces, 
respectively,  for  all  cases  of  design. 

The  units  used  for  stress  distribution  in  Figures  2.13,  2.14,  2.17,  and  2.18  are  mN/nim^.  Von  Mises  stresses  are 
represented  wherein  stresses  are  shown. 

Figure  2.13  shows  the  contact  stresses  for  the  pinion  of  the  generated  spiral  bevel  gear  drive  of  baseline  design  (GTDl) 
at  contact  point  17.  Edge  contact  has  been  found.  Design  GTD2  wherein  head-cutters  of  parabolic  profiles  have  been 
used  for  gear  generation  and  design  GTD3  wherein  head-cutters  of  parabolic  profiles  have  been  used  for  pinion 
generation  permits  the  avoidance  of  edge  contact  all  over  the  path  of  contact.  Figure  2.14  shows  the  contact  stresses  for 
the  pinion  of  the  generated  spiral  bevel  gear  drive  of  case  of  design  GTD3  at  contact  point  1 7  as  well.  Edge  contact  has 
been  avoided  indeed  and  contact  stresses  decreased  fi'om  1518  to  1175  MPa. 

Figure  2.15  shows  the  evolution  of  contact  stresses  (Fig  2.15(a))  and  bending  stresses  (Fig.  2.15(b))  for  the  pinion  of 
cases  of  design  GTDl,  GTD2,  GTD3.  Edge  contact  has  been  avoided  and  contact  stresses  decreased  for  cases  of 
design  GTD2  and  GTD3.  However,  case  of  design  GTD3  permits  the  larger  reduction  of  contact  stresses  during  the 
cycle  of  meshing 

Bending  stresses  are  almost  the  same  for  the  pinion  of  cases  of  design  GTDl,  GTD2,  and  GTD3.  Only  a  very  small 
reduction  of  bending  stresses  has  been  found  for  case  of  design  GTD3.  The  larger  bending  stresses  occurs  at  the  mean 
contact  point 

Figure  2.16  shows  the  evolution  of  contact  and  bending  stresses  for  the  gear  of  cases  of  design  GTDl,  GTD2,  and 
GTD3.  Reduction  of  contact  stresses  could  be  obtained  as  well  for  cases  of  design  GTD2  and  GTD3  with  respect  to  the 
baseline  design  GTDl  (see  Fig.  2.16(a)).  The  bending  stresses  for  the  gear  of  all  cases  of  design  are  almost  the  same. 
Larger  bending  stress  occurs  at  the  mean  contact  point  as  well. 

Figures  2.17  and  2.18  show  the  contact  stresses  for  the  pinion  and  the  gear,  respectively,  of  case  of  design  FTDl 
wherein  the  gear  is  manufactured  using  the  formate-cut  method.  Edge  contact  has  been  avoided  for  the  whole  path  of 
contact  for  one  cycle  of  meshing. 

Figure  2.19  shows  the  evolution  of  contact  and  bending  stresses  for  the  pinion  of  cases  of  design  FTDl,  FTD2,  and  the 
baseline  design  GTDl.  Contact  and  bending  stresses  for  cases  of  design  FTDl  and  FTD2  have  been  reduced  with 
respect  to  the  baseline  design  GTDl .  Also  contact  and  bending  stresses  are  lower  for  case  of  design  FTDl  with  respect 
to  the  case  of  design  FTD2  wherein  a  smaller  contact  ellipse  has  been  considered  to  increase  the  tolerances  to  errors  of 
alignment. 

Figure  2.20  shows  the  evolution  of  contact  and  bending  stresses  for  the  gear  of  cases  of  design  FTDl,  FTD2,  and 
GTDl.  Contact  stresses  obtained  for  the  gear  of  case  of  design  FTDl  are  similar  to  those  obtained  for  the  baseline 
design  GTDl.  Case  of  design  FTD2  yields  larger  contact  stresses  (see  Fig.  2.20(a)).  Bending  stresses  are  lower  for  the 
gear  of  case  of  design  FTDl . 

Based  on  the  obtained  FEA  results,  the  following  conclusions  can  be  obtained: 

Application  of  head-cutters  of  parabolic  profiles  for  generation  of  the  gear  or  the  pinion  of  face-milled  generated  spiral 
bevel  gear  drives  permits  the  avoidance  of  hidden  areas  of  edge  contact  for  high-loaded  gear  drives. 

The  application  of  head-cutters  of  parabolic  profiles  for  generation  of  the  pinion  is  more  effective  than  its  application 
for  generation  of  the  gear  for  avoidance  of  edge  contact  and  reduction  of  contact  stresses  (see  Figure  2. 1 5(a)). 

Using  the  formate-cut  method  and  a  head-cutter  of  parabolic  profile,  a  suitable  design  could  be  found  for  the  spiral 
bevel  gear  19  x  62  with  longitudinal  path  of  contact,  no  edge  contacts,  and  low  level  of  contact  and  bending  stresses. 

Designs  FTDl  and  FTD2  based  on  the  formate-cut  method  for  grinding  of  the  gear,  yield  lower  level  of  contact  and 
bending  stresses  with  respect  to  the  baseline  design  based  on  the  generated  method  for  the  gear  (see  Figures  2.19  and 
2.20). 
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Design  FTDl  wherein  a  larger  contact  ellipse  than  the  one  corresponding  to  design  FTD2  was  designed  yields  the 
lower  contact  and  bending  stresses  (see  Figures  2.19  and  2.20). 


Figure  2.13.  Contact  stresses  for  the  pinion  at  contact  point  1 7  for  baseline  design  (GTDl). 


Figure  2.14.  Contact  stresses  for  the  pinion  at  contact  point  17  for  case  of  design  GTD  wherein  the 
pinion  is  generated  with  head-cutters  of  parabolic  profiles. 
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Figure  2.15.  Evolution  of  (a)  contact  stresses  and  (b)  bending  stresses  for  the  pinion  of 
cases  of  design  GTD 1 ,  GTD2,  GTD3 . 
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Bending  Stress  (MPa)  Contact  Stress  (MPa) 


GEAR  CONTACT  STRESSES 


GTDl  -*-0702  -A-GTD3 


(a) 

GEAR  BENDING  STRESSES 


(b) 

Figure  2.16.  Evolution  of  (a)  contact  stresses  and  (b)  bending  stresses  for  the  gear 
of  cases  of  design  GTDl ,  GTD2,  and  GTD3. 
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Bending  Stress  (MPa)  Contact  Stress  (MPa) 


PINION  CONTACT  STRESSES 


-♦-FTD1  -*-FTD2  -^GTD1 


(a) 

PINION  BENDING  STRESSES 
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Figure  2.19.  Evolution  of  (a)  contact  stresses  and  (b)  bending  stresses  for  the  pinion 
of  cases  of  design  FTDl,  FTD2,  and  baseline  design  GTDl. 


NASA/CR— 2003-2 1 2336 


52 


Bending  Stress  (MPa)  Contact  Stress  (MPa) 


GEAR  CONTACT  STRESSES 


GEAR  BENDING  STRESSES 


Figure  2.20.  Evolution  of  (a)  contact  stresses  and  (b)  bending  stresses  for  the  gear 
of  cases  of  design  FTDl,  FTD2,  and  baseline  design  GTDl. 
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Chapter  3 

Conclusions 


Based  on  the  research  performed,  the  following  conclusions  can  be  obtained: 

•  A  new  computer  program  has  been  developed  that  enables  to  provide  optimized  machine-tool  settings  for 
finishing  of  spiral  bevel  gears  utilizing  the  "Formate"  grinding  process.  Such  settings  provide  minimized 
transmission  erroi^  at  a  low  level  and  a  stable  bearing  contact  adjusted  for  ^signed  misalignment.  The 
computer  program  developed  is  based  on  the  following  b^ic  ideas: 

o  Application  of  a  parabolic  fimction  of  transmission  errors  for  reduction  of  noise  and  vibration. 

o  Optimization  of  machine-tool  settings  for  the  formate  grinding  method,  including  application  of 
parabolic  cutter  blade  profiles. 

o  Application  of  an  enhanced  approach  for  automatic  development  of  finite  element  models  and  stress 
analysis 

•  The  developed  approach  of  design,  generation,  simulation  of  meshing  and  contact,  and  stress  analysis  of 
formate  cut  spiral  bevel  gear  drives  has  been  successfully  applied  for  the  design  of  a  formate  cut  spiral  bevel 
gear  drive  with  a  stabilized  bearing  contact  and  favorable  shape  of  transmission  errors  of  low  magnitude. 

•  A  formate  cut  spiral  bevel  gear  drive  has  been  designed  considering  as  a  reference  a  design  b^ed  on  face- 
milled  generated  spiral  bevel  gear  drives.  The  provided  design  allows  to  increase  the  productivity  of 
manufacturing  of  the  gears  by  the  formate  cut  method,  avoid  edge  contacts,  and  reduce  contact  and  bending 
stresses. 

•  Head-cutters  with  parabolic  profiles  have  been  applied  for  the  generation  of  the  gear  of  formate  cut  gear 
drives  and  for  generation  of  pinion  and  gear  of  face-milled  generated  spiral  bevel  gear  drives.  Application  of 
head-cutters  with  parabolic  profiles  allows  us  to  increase  the  mismatch  between  generating  surfaces  of  pinion 
and  gear  and  avoid  the  appearance  of  edge  contact. 

•  Automatic  generation  of  finite  element  models  have  been  developed  to  perform  multi-tooth  stress  analysis  for 
several  contact  points  along  the  path  of  contact  and  obtain  the  evolution  of  contact  and  bending  stresses. 
Hidden  areas  of  edge-contact  for  high-loaded  gear  drives  are  detected  if  they  exist. 

•  A  method  for  compensation  of  the  shift  of  the  path  of  contact  caused  by  a  change  of  the  shortest  distance 

between  the  axis  of  the  pinion  and  the  gear,  AE ,  has  been  proposed.  It  is  b^ed  on  the  compensation  of  error 
of  alignment  AE  by  the  axial  displacement  of  the  pinion,  . 
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